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Abstract 

We consider asymptotics of the correlation functions of characteristic polynomials 
corresponding to random weighted G(n, Erdos - Renyi graphs with Gaussian weights 
in the case of finite p and also when p —>■ oo. It is shown that for finite p the second 
correlation function demonstrates a kind of transition: when p < 2 it factorizes in the 
limit n oo, while for p > 2 there appears an interval (—A*(p), A,(p)) such that for 
Ao G (—A*(p), A*(p)) the second correlation function behaves like that for GUE, while 
for Ao outside the interval the second correlation function is still factorized. For p —>■ oo 
there is also a threshold in the behavior of the second correlation function near Aq = ±2: 
for p <C the second correlation function factorizes, whereas for p ^ it behaves 
like that for GUE. For any rate of p —oo the asymptotics of correlation functions of any 
even order for Aq £ (—2, 2) coincide with that for GUE. 


1 Introduction and main results 


Consider an ensemble of hermitian n x n random matrices of the form 




( 1 . 1 ) 


where 


djk = P 


1 with probability 


P 


n 


0 with probability 1 — 


P 


n 


(1) , • (2) • / 7 


( 1 . 2 ) 


and {wp} ,wii :l<j<k<n,l<l<n} are i.i.d. random variables with zero mean 
such that 

2E{|u;W| 2} = 2E{|u;jyId = B{\wa\^} = 1, j y k. (1.3) 

Here and everywhere below E denotes the expectation with respect to all random variables. 
{djk ■ j < k} are also independent of each other and of ,wii. 

These matrices are known as “weighted” adjacency matrices of random Erdos - Renyi 
G(n,^) graphs, with {djk} corresponding to the standard adjacency matrix and {wjk} — 
the set of independent weights, which we take to be Gaussian. These matrices are widely 
discussed in the last few years since they demonstrate a kind of interpolation between a 
“sparse” matrix with finite p, when there is only a finite number of nonzero elements in each 
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line, and the matrix with p = n coinciding with Gaussian Unitary Ensemble (GUE). The 
results on the convergence of normalized eigenvalue counting measure 

Nn{A) = #{xf G A, j = 1,..., n}/n, iV„(R) = 1 

of these matrices in the case of finite p were obtain in [201 |2T] on the physical level of rigour, 
then in [T] for Wjk = 1 and in |14j for arbitrary {wjk} independent on {djk} and having 4 
moments. 

It was also shown that for p —oo the limiting eigenvalue distribution coincides with GUE 
Urn Nn{A) = [ psc{X)dX, PscW = • 1[_2,2], (1-4) 

n—>oo J ZTT 

A 

while for finite p the limiting measure is a solution of a rather complicated nonlinear integral 

equation which is difficult for the analysis. It is known that the support of the limiting 

measure (spectrum) is the whole real line. Gentral Limit Theorem for linear eigenvalue 

statistics was proven in [26] for finite p and in m for p —>• oo. Eor the local regime it 

was conjectured the existence of the critical value pc > 1 (see 0) such that for p > Pc the 

eigenvalues are strongly correlated and are characterized by GUE matrix statistics, for p < Pc 

the eigenvalues are uncorrelated and follow Poison statistics. The conjecture was confirmed 

by numerical calculations m and by supersymmetry approach (SUSY) [8l|T8| on the physical 

lever of rigour. Notice, that the results of the present paper confirm the existence of similar 

threshold for the second correlation function of characteristic polynomials. Rigorous results 

for the local eigenvalue statistics were obtained recently in lattoj. First for p ^ and 

then for p^ rf with any e > 0 it was shown that the spectral correlation functions of sparse 

hermitian random matrices in the bulk of the spectrum converge in the weak sense to that 

of GUE. For the edge of the spectrum, it was proved in [T2| that for p S> the limiting 

probability PlmaxA^”^ > 2 + admits a certain universal upper bound, whereas the 

j ^ 

result of m implies that for p <C the limiting probability P{maxA^"'^ > 2 + x/p} is zero. 

Note that more advanced results for the edge eigenvalue statistics were obtained in [28| for 
so-called random d-regular graphs. It was shown that if 3 < d <C and Wjk = ±1 then 
the scaled largest eigenvalue or (HI) converges in distribution to the Tracy-Widom law. 

The correlation functions of characteristic polynomials formally do not characterize the 
local eigenvalue statistics. However, from the SUSY point of view, their analysis is similar 
to that for spectral correlation functions. In combination with the fact that the analysis for 
correlation functions of characteristic polynomials usually is simpler than that for spectral 
correlation functions, it causes that such an analysis is often the first step in studies of local 
regimes. 

The moments of the characteristic polynomials were studied for a lot of random matrix 
ensembles: for Gaussian Orthogonal Ensemble in |6|, for Gircular Unitary Ensemble in ffUE], 
for /3-models with /3 = 2 in [2l|29] and with /3 = 1, 2, 4 in [3l[l7]. Gotze and Kosters in |9| have 
studied the second order correlation function of the characteristic polynomials of the Wigner 
matrix with arbitrary distributed entries, possessing the forth moments, by the method of 
generation functions. The result was generalized soon on the correlation functions of any 
even order by T. Shcherbina in |22| where it was proposed the method which allowed to 
apply SUSY technique (or the Grassmann integration technique) to study the correlation 
functions of characteristic polynomials of random matrices with non Gaussian entries. The 
proposed method appeared to be rather powerful and since that was successfully applied to 
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study characteristic polynomials of sample covariance matrices (see [23]) and band matrices 

(|M1[25|). 

In the present paper we apply the method of m to study characteristic polynomials 
of sparse matrices. To be more precise, let us introduce our main definitions. The mixed 
moments or the correlation functions of the characteristic polynomials have the form 


2m 


F2m{A) = 'E\ lldet{Mn-Xj 


i=i 


(1.5) 


where A = diagjAi,..., A 2 m} are real or complex parameters which may depend on n. 

We are interested in the asymptotic behavior of (II.51) for matrices (|l.ip . as n —oo, for 

^ • _ 

Aj = Ao + —, j = l,2m, 
n 

where Aq, are real numbers and notation j = 1,2m means that j varies from 1 to 

2m. 

Set also 


T>2m(A) 


F2m{A) 



A*(j?) 


\/4-8/p, if p > 2; 
0, ifp<2. 


( 1 . 6 ) 


Theorem 1 Let an ensemble of sparse random matrices he defined by (|l.ll) - (|1.3p for finite p 
and let be Gaussian random variables. Then the correlation function of two char¬ 

acteristic polynomials (lESl) for m = 1 satisfies the asymptotic relations 

(i) for Ao E (-A*(p), A*(p)) 


lim D 2 (A) 

n^oo 


sin((xi - X 2 )^/A^(p)^ - A^/2) ^ 
(xi - X2 )\/A*(p)2 - A^/2 


(ii) for Ao ^ (-A*(p), A*(p)) 


lim D 2 (A) = 1, 


where D 2 and A* (p) are defined in ()1.6p . 

Remarks 

1. The theorem shows that the second order correlation function has a threshold p = 2, 
i.e. if p > 2 there are two types of the asymptotic behavior — cases (i) and (ii), if p < 2 
there is only one type of the asymptotic behavior — case (ii). 

2. If we let Ao depend on n, the asymptotic regimes (i) and (ii) are fully agreed. 

3. Note that A*(p) —>■ 2, as p —)• 00 , and since the limiting spectrum is always [—2, 2] (see 
(El), therefore for p ^ 00 one expects GUE behavior for all Ao E (~2, 2) (i.e. for all 
Aq in the interior of the limiting spectrum, cf (ll.4p i: we confirm this in Theorem [2j 
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Theorem 2 Let an ensemble of diluted random matrices he defined by (jl.lh - (ll.3p . p —>• oo 
and let be Gaussian random variables. Then the correlation function of character¬ 

istic polynomials (113]) for Ao G (—2, 2) satisfies the asymptotic relations 


lim -D 2 m(A) 

n—>oo 


S2m{X) 
S2m{I) ’ 


with X = diagjxi,..., X 2 m} cind 


det 


S2m{X) = 


sin(7rp^e(Ao)(Xj - Xm+k)) y^ 
^Psc('^o) ^m+k) f j fc=l 


A(xi, . . . , • • • ) X2m) 

where A(yi,..., ym) is the Vandermonde determinant of yi,... ,ym. 


(1.7) 


Notice that S 2 m{I) is well defined because the difference of the rows ji and j 2 in the 
determinant in (II.7p is of order — xj^), as Xjj —)• xj^- The same is true for columns. 

To formulate our last result we introduce the Airy kernel 

x-y 

where Ai{x) is the Airy function. 


Theorem 3 Let an ensemble of diluted random matrices he defined by (ILI|)-(II3|), p 
oo, and let Gaussian random variables. Then the correlation function of two 

characteristic polynomials (HSI) for m = 1 satisfies the asymptotic relations 


(^) 




lim Zl 2 ( 2 / + n"^/^A) = 1; 

n—¥oo 


(it) If^^c 

lim D 2 { 2 I + n-‘^/^X) 

n^oo 


A(xi + 2c, X 2 + 2c) 

■^A(xi + 2c, xi + 2 c)A(x 2 + 2c, X 2 + 2c) ’ 


where D 2 is defined in (HH) and A is defined in dEll). For Ao = — 2 similar assertions are 
also valid. 


Remarks 

1. Notice, that the case p S> corresponds to the case c = 0 in (ii). 

2. The results of Theorem [3| are in a good agreement with the results of [121113] in the 

sense that the asymptotic behavior changes when p crosses the rate However, in 

m it is argued that in the case p <C A the appropriate scale is p ^ instead of n . 
We postpone the study of F 2 with the scaling p~^, as well as the related study of F 2 
near A*(p) for finite p, to subsequent publications. 

The paper is organized as follows. In Section [2] we obtain a convenient integral representa¬ 
tion for F 2 m using integration over the Grassmann variables and Harish-Chandra/Itzykson- 
Zuber formula for integrals over the unitary group. Sections [3] jj] and [5] deal with the proof 
of the Theorems [U [2] and [3] respectively. The proof is based on the steepest descent method 
applied to the integral representation. 

Notice also that everywhere below we denote by C various n-independent constants, which 
can be different in different formulas. 
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2 Integral representation 


To formulate the result of the section, the following notations are introduced 

• A(diag{yj}^^^) = ^{{yj}j=i) is the Vandermonde determinant of {yj}j^i] (2-1) 

• End V is the set of linear operators on a linear space V ; 

• = {a G < ai < ... < afc < n}. (2.2) 

The lexicographical order on In,k is denoted by 

• is the space of self-adjoint operators in EndA^C^™' (see [301 Chapter 8.4] for 
definition of A'^V); 

• dB = dPi{B) = n dBaa n dSiBapd^Bap (2.3) 

a&hm.l OL^P 

is a measure on Bap denotes the corresponding entry of the matrix of B in some 

basis. It is easy to see that dPi{UBU*) = dPi{B) for any unitary matrix U, so the 
definition is correct. 

2m 

• H^=Y\ (2.4) 

1=2 

• Set also 


2m 


2m 


fcl+2/c2+---+2mfc2m=2m ^ ^ 


Aks 


(2.5) 


where are the sequences of certain n,p-dependent numbers and 

G = {G2,...,G2m), Gi GEndA'C^"^, f 

Exterior product A A B of operators is defined in Section [6.11 Since dimA^”*C^”^ = 1, 

_ 2m 

the space EndA^’^C^’” may be identified with the C. In (12.5p /\ (bgGs — is 


S = 1 


understood as s /\ [bgGg — bgl)^^” ^ 

• = TT™ exp g x]]. 

'' j=l ■> 

We prove the following integral representation for the correlation function F 2 m- 


( 2 . 6 ) 


Proposition 1 Let Mn be a random matrix of the form (|l.ip - (|1.3l) . where j < k, 

wii have Gaussian distribution. Then the correlation function (HSI) admits the representation 


j2m? m ex;p 


F2™(A) = 0(2™) (V)- 


^ 2m N 

{^0 E Xj\ 
r j=l > 


A{X) 


Hm R2m 


^ 2m N 

A(r)expi -iJ2xjtj\e^B^^BR)dTdR, (2.7) 
^ i=i ^ 
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2m 

where R = {R 2 ,, R 2 m), Ri G R 2 m,i, dR = Y\ dRj, T 

i=2 


diag{i,}2-, dT 


2m 


n dtj, 

i=i 


/ 2 m(T, R) = log A 2 m{T, R) 



and all other notation is defined at the beginning of the section. 


( 2 . 8 ) 


Remark 1 In the special ease m = 1, the representation (1^ simplifies to 


F2{K)=Cn{X) 


where Cn{X) = Cn\x) and 


le 


Ao {X1+X2) 


Xi - X2 


/ ih - t 2 ) exp \ Xjt 

i ^ i=i 


(2.9) 


f{T, s) = log{b2S - tit2) - ^ + *^o)^ + 


( 2 . 10 ) 


The proof of Proposition [T] is based on the method of integration over the Grassmann 
variables, the required properties of which are reviewed in Section 16.11 The proof of the 
proposition is given in Section 12.11 


2.1 Proof of Proposition [T] 

Let us transform F 2 m{X), using (16.ip 


f 2m .. ( r ( / 

T 2 m(A) = pi det(M - Aj) I = / exp I - ^ f ^ djkWjkfi>jif^ki 

^ j=l ^ ^ d I V i<j^k<n 




2m n 


Y. u n n dipjidipji 


1 = 1 j=i 


Averaging first with respect to {wjk}, we obtain 


1 


2m 


2m n 


F 2 m{X) — E< / exp < ^ djkXjkXkj +E U44+E nn difjidfiji L 

^ j<k j=l ^ 1=1 ^ ^ «=li=l ^ 

where, in order to simplify formulas below we denote 


2m 

Ajfc = 

1=1 


Since evidently {XjkXkjY'^~^^ = Oj we get 

✓ 2m . X 2m .. .. 

E } = Ei 1 + J]; dfkixjkXkjY [ = 1 + Z] ^ • FTFfii^FXkjy, j < k. 

1=1 ^ i=i ' ^ 
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Define the numbers by the identity 


y- 2m N 2m ^ ^ 

exp '^aiixjkXkjY f = 1 + yy • ZJZTziXjkXkjY- 

^ 1=1 J 1=1 p ^ 


Observe that 


and that 


1 n — p 

oi = —, 02 = —T 

n Zpn^ 


ai 


C 


pi 


n ^ oo, 


ai = 0, I > 1, if p = n. 


Then F 2 m(A) can be represented as 

2m 


/ y- 2m n / ^ \ 

exp EE ^liXjkXkj) 

^ jck l=l j=l ^ l=l ^ 

2m n 

nn diijidijjji. 

1=1j=i 

To facilitate the reading, the remaining steps are first explained in the simpler case 
and only then in the general case. 


2.1.1 Case m = 1 

Let us transform the exponent of (|2.12p . 

XjkXkj = '^jai'^jf}i'>PkPi'>Pko 

“i/3S/2,1 0,/3S/2,1 

2 2 2 2 

iXjkXkjf = ‘iYl'lpjli^kl'tPkl'^jl = 4:Yl'tpjliljl'tpkl'4^kl = ‘iYl'^Pjl'tl^jlYl'lPkq'^Pkq, 

1=1 q=l 


kai ; 


1=1 


1=1 


where / 2 ,i is defined in (12.2D . Since = pip = 0, we have 


XjkXkj + Y = -Y\{Y • 

j<k j=i 1=1 ^j=i ^ ^ j ^ ^ j ^ 


■ ]=1 3 

n 2 \ 2 


j<k ^j=l1=1 

Hubbard-Stratonovich transformation (16.2D applied to (I2.12D yields 


exp <! oi ( XjkXkj + Y I = ” 


j<k 


i=i 


2^^V 27r2 


/ 

W2 


expi - - I + ' 

i=i 


n exp ^ Z] ti'ipji'ipji + i{u - iv)'ipj2'4’ji + iiu + iv)'tpj^-ipj2 \dQ; 

j=i ^ 1=1 


( 2 . 11 ) 


( 2 . 12 ) 
m = 1 
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X n exp { -s^/4na2 ■ il^ji'ipji'ipj2'ipj2} ds, 
i=i 

where "^2 is the space of self-adjoint operators in EndC^ and 

n — ( tt -I- if V 

\ u-iv t 2 y ’ 

dQ = dtidt2dudv. 


Set 62 = —^/^na 2 = — . Now we can integrate over the Grassmann variables 

V 

e-t"' j{b 2 S - det {Q - iA))"e-t t'-Q'dQds 

R H 2 

Change the variables tj —>■ tj + iXj and move the line of integration back to the real 
axis. Indeed, consider the rectangular contour with vertices in the points (—i?, 0), {R,0), 
{R, —iXj) and {—R, —iXj). Since the integrand is a holomorphic on C function, the integral 
over this contour is zero. Because the integrand is a polynomial multiplied by exponent, the 
integral over the vertical sides of the contour tends to 0 when i? —)• 00 . So, recalling that 
Xj = Xq + Xj/n, we can write 




F 2 { A ) = . £-^0(2:1+2:2) 


vr 


X / e 2^ 


' J ( 62 s - detQ)’"e-tt''exp{-itr XQjdQds, 
H 2 


where 


\ 3/2 

a(7 = n(|-J 


(2.13) 


Let us change the variables Q —>■ U*TU, where U is a unitary matrix and T = diagjti, t 2 }- 
Then dQ changes to 5(ti — t 2 )‘^dtidt 2 and 


F2{A) 




J Jih-t2)\b2S-tR2r 

R 

U 2 




where U 2 is the group of the unitary 2x2 matrices, dU 2 {U) is the normalized to unity Haar 
measure, dT = dtidt 2 - 

The integration over the unitary group using the Harish-Chandra/Itsykson-Zuber formula 
(j6.3l) implies the assertion of Remark [TJ 
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2.1.2 General case m > 1 

Let us transform the exponent of (|2.12p . 


iXjkXkj) — (^0 E n '^jag'^kaq^’kPq'^jPq 

a,P&l2m,l 9=1 

I I 

0./3G/2m,i 9=1 


9=1 


where l 2 m,i is defined in (|2.2I) . 


_ 2 

Since V’|; = ipji = 0, we have 
iXjkXkj) + 

j<k 


Y^^xjkXkjy+ ^ nn 

i=l ^i<fca,/3e/2m.i 9=1 9=1 


1 


n I 


I 


+5En '^jaq 

i=l 9=1 9=1 


I I 

(-i)'(i!)" E ED i^jaq n '^kPq'^kaq 

«,/ 3 e/ 2 m,i j,k 9=1 9=1 

(-i)'(i!fE(Eri ' 4 ’jaq )(En '^kPq'll^kaq 

a,13 ^ j q=l ^ k q=l 

n I 




I^jaq 


a ''j=lq=l 

( n I \ y n I 

E ) (E ri^i/3.^. 

uL-,^ j = l q=l 


j=lq=l 


9«9 / ’ 


where -< is the lexicographical order. 

Hubbard-Stratonovich transformation (|6.2p yields 


exp <1 Yl 


n I 


JO^q 


j = l q=l 



^ / exp| 6 ;(Qi)«af j “ 

'9=19=1 


n I 


n 


n I 


exp](-l)'(«fa,( 5 ;n ‘^joiq'^jPq En 

■ 9 = 19=1 ' ^ 9 = 19=1 


9«? 


re 

vr 


exp{-re|(Qz)„/3p} i ( E 11 ) (Q0«/3 

^i=lg=l ^ 


+ &Z (En ' 4 ^jfiq'^joiq j (Qi)a /3 f )a/ 3 ) 

■ 9 = 1 9=1 
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where 


bi = iHly/imi. (2-14) 

The above computations compose into the following representation of the exponent of ()2.12p 


exp i ai ( Y^iXjkXkjy + 


j<k 


i=i 


l/ 2 m\ 

2 \ I ) / fi\ 2 \ I ) 
TT 


X [ exp|-^trQf| Jlexpjft/ 

nt ^ ^ 

where 'H 2 m,i is the space of self-adjoint operators in EndA^C^™ and dQi is defined in (12.3p . 
Therefore, substitution of (I2.15P into (|2.12l) gives us 

/ 2m n f I 

n qf } n n 


where is defined in (I2.4D and 


1 1 
2 




2m 


^ 1=1 




TT 


Now we can expand the exponents of (j2.16p into the series 

2m I 


1 1 


exp 

1=1 a,j3 q=l 

= exp|^^(4)„/3n^ja,V'i/3d ’ (2.17) 


1=1 a,, 


q=l 


2m ^ / 2m 


k\ , 

k=l ^ 1=1 q:,/3 


q=l 


where _ ^ ^ 

Qi = biQi + A, Qi = biQi — bil, 

hi = ( 2 “^ — 2~’')ai, hi-i = 0 . 

The most important terms contain all 4m Grassmann variables {ipjs, V'jsiF™!; because the 
other terms become zeros after integration over Grassmann variables. Thus, expansion of 
p2.17p with Lemma [6] implies 


2m 


®^p \ Yi n \ n 


1 = 1 Oi,^ 


q=l 


2m 


1=1 


E 


{ki -|- ... -|- A: 2 m)! 


(2m) 




kj£Z+ 


X E f A ^ 

G()/3G/2m,2m ^—1 


2m \ 2m 2m 

q=l 1=1 
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where only the most important terms remain. Integration over the Grassmann variables and 
substitution of the result into (|2.16p gives us 


i"2m(A) 


4'™^/ j {A2m{Qi,Q)y 

HmH2 



n 

2 


2m V 

^trQf ^QidQ, 


1=1 


y. 2m 

where Qi = Qi + = Qi - iA, Q = {Q 2 ,, Q 2 m), dQ = Y\ dQi and ^ 2 m is dehned in 

1=2 

dZS]). 

Change the variables {Qi)jj {Qi)jj + and move the line of integration back to the 
real axis. Similarly to the case m = 1, the Cauchy theorem yields 


A2m(A) 




vr" 


• exp 


2m V „ 

7 = 1 ^ . 


{A2m{Qi,Q)y 


^ Hm W2m,l 

exp I - ^ ( tr(Qi + iAo)^ + ^ tr Q' 


2m 


1=2 


^—itrXQi 


dQ 1 dQ , 


where Cn'^\x) is defined in (j2.6[) . Let us change the variables Qi = U*TU, where U is a 
unitary operator and T = diag{fj}|™^. Then dQ changes to 7r”^iL^A(T)^(ir and 


F2m{A) — 


d^2m 


exp{A„£x,}/ / 

3 1 ^ Tm9m 


X{Tf {A2m{QuQ)y 


I' / 2m 2m \ ^ „ 

exp{+ J 

^ = ^ ' = 2 U2m 


-itrXU*TU 


dU2m{U)dTdQ, (2.18) 


2m 

where K 2 m = n j!, A(T) is dehned in dZII), U 2 m is the subgroup of the unitary operators 
i=i 

2m 

in EndC^”^, dU 2 m{U) is the normalized to unity Haar measure, dT = ])([ dtj. Transform 

i=i 

A2m(<3lj Q) 

2m 2m 

{biQiT^^ A l\{bsQs - bsl)^^‘ = {WbiTUf'^^ A l\{{U*Uf\hsQs - bsI){U*U)^T^‘, 

s=2 s=2 

where I is the identity operator. The assertion (iv) of Proposition [2] implies 


2m 


\As\Akis 


TA2m 


{U*biTU)^'^^ A /\{{U*U)^^{bsQs -bsI){U*Uy 

s=2 

z 2m 

= ([/*)^ 2 m J ^ (bsU^VsiU*)^" - bsl)^^^ ] 

2m 

= {biT)^’^^ A /\{bsU^^QsiU*)^^ - bslT^‘. 


s=2 


s=2 
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Change the variables Qi = RiU^\ I = 2 , 2m. Since 17'^^ is the unitary operator, dQ 

2m 

changes to dR = H dRi. Then (I2.18P implies 
1=2 


T2m(A) = 


dn'"\X) 

I^2m 


exp{A„£x,} / / 

.7 — 1 U TrB9m 


AiT)^{A2miT,R)y 


3 = 
2m 


Hm 


f / zrn 2m \ x « 

xexp{-^(^^(tj+ao)2 + ^tri?fj} J 

where R = {R 2 , ■ ■ ■, R 2 m)- 

I iA2miT,R)) 


U2r, 


-itTXU*TU 


dU2miU)dTdR, 


exp 


H„ 


ti' Rf ] ^dR 

j=i 1=2 


(2.19) 


is a symmetric function of Indeed, after swapping and and changing the 

variables Ri —)• where is the unit matrix, in which rows ji and j 2 

are swapped, the integrand in (|2.19p remains unchanged. Hence, the integration over the 
unitary group using the Harish-Chandra/Itsykson-Zuber formula (j6.3p can be done, which 
yields the assertion of Proposition [TJ 


3 Proof of Theorem [T] 


To find asymptotics of T 2 (A), we apply the steepest descent method to the integral represen¬ 
tation As usual, the key technical point of the steepest descent method is to choose a 

good contour of integration (in our case it is 3 dimension space of {ti,t 2 ,s)), which contains 
the stationary point s*) of / then to prove that for any (ti,t 2 , s) in our “contour” 


^f(tut2,s)<9tf(tl,tys*). 


(3.1) 


Let us introduce the function hr 


1 / ^ 

ha(tl,t2,S,b2,Ao) = 9 ( logA - 


2 2 

62 - 2 a(l - q;)Ao ), 


i=i 


a 


where 


A = ( 62 'S — tit2 + ct^Aq) + o3 X^iti + t 2 )^- 

Then ??/ at our “contour” (which is defined further) has the form 

3?/(r,s) = h„(5Rfi,m2,s,62,Ao) + - bl + il-a)Xl 

for some a. To prove (j3.ip . we use the following lemma. 


(3.2) 

(3.3) 


Lemma 1 Let ha be defined by (13.2p and (13.3p . Then for every a G [1/2, 1), ti, t 2 , s, b 2 , 
Ao G K the following inequality holds 

h«(ti,t 2 ,s, 62 ,Ao) < ^log -1 (3.4) 

Moreover, the equality holds if and only if at least one of the following conditions is satisfied 
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(a) a = 1 / 2 , ti = —t 2 = ±a/4 — 463 — Aq/ 2 , s = 62 ; 

(b) a = 1/2, ti = t2 = ±^4 - 46 ^ - Aq/2, s = -62, &2A0 = 0; 

/c/ ti = t 2 = 0, s = 62 ^, a(l - a)A§ + = 1 ; 

(d) ti =t 2 = 0, s = -b 2 ^, 62 = ± 1 ^, Ao = 0 . 

Proof. Rewrite the inequality (13.41) in the form 


log 


1 — a 


a 


^ 1/2 < 1 ^^2 ^2 g 2 ^2 2a(l - a)Ao) , 


where d = ^—^ 62 . Since 

rv " 


logL^AV2<L^^l/2_i^ 


a 


a 


it is sufficient to prove 

A <^(tl + tl + + d"^ + 2a(l - a)Ao)^ • 

Recalling (j3.3p . we have 


(3.5) 


(3.6) 


1 — a 


a 


A = s^d^ + 


1 — a 


a 


tft2 + a^(l — a)^AQ 
1 — a 


- 2 - 


a 


sdtit2 + 2a{l — a)XQsd + (1 — Q:)"^AQ(q — t^). 


dM]) is transformed into 
2 


s^d^ + 


1 — a 


a 


tfti - 2- 


1 — a 


a 


sdtit2 + 2 a(l — a)XQsd 


+ (1 - a)2A§(tf - ti) < l(t; +ti)2 + f 

+ 2 (^1 ^2)('®^ ffi ^^) ffi ~ '^)Ao(^i ffi ^ 2 ) ffi ~ '^)Ao(®^ ffi d^). 


The last inequality is the sum of following obvious inequalities 
(1 — a)^Ao(ti — tl) ^ 0^(1 “ ®)Ao(^i + ^ 2 )) 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

2 a(l — a)Aosd < a(l — a)Ao(s^ + df). (3-11) 

It remains to determine conditions when the equality in (13.41) holds. It holds if and only 
if the equalities in (|3.5I) . (I3.7l) - (|3.1ip hold. Let (n') denotes the corresponding equality for 
inequality (n). Then 


\ — a 
a 

1 — a 


s 2 ^ 2 < 1(^2 ^^ 2 ) 2 ^ 

1 


-2- sdtit2 < 2\sdtit2\ < -(tf + + d ), 

a 2 

2„j ^ „ /I \2/„2 , i2\ 


(1331 1 4^ A^/"^ = 


a 


1 — a ’ 


(3.12) 


^ = df, 

^tl= tl 
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Everywhere below until the end of the proof we assume that = d? and t\ = t^- Then 


(1321) ^ 


dSIB) ^ 


a = 1/2; 

Aoti = 0; 
a = 1 / 2 ; 
ti = 0; 

fa = 1/2; 

< 0 ; 

tis = 0 ; 

sd > 0 ; 

Ao = 0. 


Let us consider the following cases 

1 . h = 0 . 

1.1. Ao = 0. Then ()3.12p is transformed into {b 2 s)‘^ = Since we get 

62 = that implies (d). 

1.2. sd > 0. Then (I3.12p is equivalent to a^Ag + ^-^&2 = that implies (c). 

2 . 7 ^ 0 => a = 1 / 2 . 

2.1. s = 0. Hence, (|3.12p is transformed into Aq/4 + t\ = 1 that implies (b). 

2 . 2 . s 7 ^ 0 ^ dstit 2 < 0 . 

2.2.1. sd > 0. Then ()3.12l) is transformed into Aq/4 + 62 + = 1- Condition (a) is 

satisfied. 

2.2.2. sd < 0 ^ Ao = 0. Then (j3.12p is transformed into = 1. Condition (b) 

is satished. 


Finally, it is easy to check that the values of ha at the points satisfying (a)-(d) are equal 
to the r.h.s. of p.4p . ■ 

Now we are ready to prove Theorem [H We start from the lemma 


Lemma 2 Let all conditions of Theorem[I\ are hold and Ao G (A*(p), A*(p)). Then T 2 (A) 
satisfies the asymptotic relation 


F 2 {A) = 2nexp{n(Ao + 


2)/2 + Ao(xi+X2)/2} 


X 


sin((a;i -X2)V4-46| - A^/2) ^ ^ 

(xi - X2) 


(3.13) 


where 62 is defined in (12.1411 . 
Proof. Set 


= diag {t*,tl} - iAo/ 2 , 


(3.14) 


where r^, = 1 , 2 . 
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Consider the contour = —Ao/2, s G M. It contains the points — 1 X 0 / 2 , — 

iAo/2,62) and {t *2 — i\o/2,t\ — iAo/2,62), which are the stationary points of /. The contour 
may contain another stationary points of /, but this fact does not affect the proof, except the 
case Ao = 0 for which the points 62) and { 1 * 2 , 12 ,—b 2 ) are also under consideration. 

First, consider the case Aq 7^ 0. 

Shift the variables trj ^ — iXo/2 in (|2.9p and restrict the integration domain by 

Br = {{ti,t 2 ,s) G : max{|ti|, |t2|, |s|} < -R} • 

Then 

Ap\o(xi+X 2)/2 f r i \ o 

F 2 {A) = Cn{X)— -— / g{T, s)e^B'^-2^°^^)dTds + 0{e-^^ /^), n ^ 00, (3.15) 

(Xi - X2) J 
Br 


where / is defined in (12.101) and 


g{T, s) = (ti - t2) exp j - * ^ Xr^t^ ^. 

77=1 

Then it is easy to see, that for g ^ R 

f{T^^^\b2) = ^bl + ^Xl-l-, 

/ (t* - iAo/2)2 + 1 bl -b2{tl-iXo/2) 

/"(ri’?"),62) = - bl {%- i\o/2f + 1 -b 2 {t* -i\o/2) 

\ -b2{tl-iXo/2) -b2{t;-iXo/2) 52 + 1 


det 

M) = 

_(4- 

-462 - Ao) < -{K{pf - 

Ag); 




/ 

K? - Ag/4 + 1 


bl 

-b2tl \ 


M) = 

- 


{t*,? - 

■A2/4 + 1 

-b2t* 



V 

- 62 C 


■b2t* 

62 + 1 / 

det3ft/"(ri^''^ 

M) = 

-(1 

-A2/4) (4-462. 

-A^X 

-{Kipf - 

- Ag)V4, 


where is defined in (j3.14p . 

Note that 

3^/ ^ ~ hi/ 2 {ti,t 2 , s, 62, Ao) + -62 + 2-^05 

where ha is defined in (|3.2I) . According to Lemma (H 3?/ (T — |Ao, s), as a function of real 
variables ti, t 2 , s, attains its maximum at {T^^''\b 2 ). Hence, 62) is nonpositive, 

but since det 3?/"62) < 0, 3?/"(Ti^^\ 62) is negative dehnite. 

Let be a Jogrr-neighborhood of the point {T^^''\b 2 ) and let Vn everywhere 

below denote the union of such neighborhoods of the stationary points under consideration, 
unless otherwise stated. Then for (T, s) ^ Vn and sufficiently large n we have 



5R/(Ti""\62)-5R/ (t-Uo,s^ 

>min min |3?/(ri’'"\ 62 ) - 3?/(t - ^Aq, s 

{T-±Ao,s)edv/’'>''^ I V 2 


>cl55T, 

n 
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Thus we can restrict the integration domain by Vn- 

Set q = q* = Then expanding / and g by the Taylor formula and 

changing the variables q n~^^‘^q + q*, we get 


FsfA) = n-^/^CJX) ■ ^ + ^2 - 2)/2 + Ao(xi + X2)/2} 

^ (xi- X 2 ) 

X J 5r(3?ri''''\62)exp|ig/"(ri’'‘'\62)g'^|dg + o(l)^. 

log n,log n]3 

This is true because g{^Tj^’^\b 2 ) / 0. Performing the Gaussian integration, we obtain 

i?2(A) = ( —(X) • ^exp{n(Ag + fej - 2)/2 + Ao(xi + X2)/2} 

\n ) {xi- X 2 ) 

X (Y,9{^T^''"\b2)det-^/\-f\T^^''\b2)} + o{l)^^ (3.16) 

Since 

62 ) = r? / I., 

and Cn has the form (I2.13p . we have (j3.13p . 

If Ao = 0, then repeating the above steps, we obtain the formula similar to (|3.16p . The 
only difference is that there are two more terms (i.e. there are as many terms as stationary 
points) in the sum. Since g is zero at the points with ti = t 2 , we have exactly (I3.16P and 
hence the asymptotic equality ()3.13p is also valid. ■ 

The assertion (i) of the theorem follows immediately from Lemma [2j 


Lemma 3 Let all conditions of Theorem{J\ are hold and Ag > 4 — + e for some e > 0. 

Then F2(A) satisfies the asymptotie relation 

(i) for Ao / 0 


T2(A) 


a^expIreAT (1-a)Ao(xi+X 2 )} . . .. 

(2a - 1)3/2(2 - a(l - a)(3 - 2a)A2)V2 


where a and A satisfy 


a G (1/2,1), 


a(l — a)Ao + b 2 — 1 = 0, 

A = f ^-zaAg, —^^2^ • 


(3.17) 


(3.18) 

(3.19) 


(a) for Ao = 0 

F2(A) = (fc + 1 + (-l)“(fe - 1)) (1 + 0(1)). (3.20) 


Proof. Choose Ati = Qt 2 = — oAq, s G M as the good contour with the stationary point 
(—mAo, —iaXo, ^^b 2 ), where a satishes (13.181) . Existence and uniqueness of such a follow 
from the fact that the l.h.s. of ()3.18p is a monotone decreasing function of a whose values at 
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a = 1/2 and a = 1 have different signs. Everywhere below we assume that a is a solution of 

(ICT) . 

If Ao = 0, we have two stationary points at the contour — (0,0, ±1). 

Consider the case Aq > e. Shifting the variables tj —>■ tj — iaXo, similarly to (13.151) we get 


It is 


„-p(l-o)Ao(a:i+a:2) f 

F2(A) = Cn{X) -^^- / g{T, /^). (3.21) 

(Xi - X2} J 

easy to check that 




/^ - iaAo, ^-^62^ = ^ + (1 - a)Ao + log 

/ 1 - 


/"(^-iaAo,t^62) 


(1^)352 i_n-n 




- a 


(3.22) 


„ , bl ib2Xo^^-^ ^ 

1 - (1 - afXl ib2Xo^^^^ 


\ zUq, zUq 2, \ oi 

det/"^ - iaAo, “^^2^ = “ a)(2Q: - l)Ao + 26| ^ 

iaAo, ^ ^ ^ ^ ^ a ) 

/ /I \2 \ 

. r./l— r>\\ 


^ + i>uXf ) 

V O V 

<0; (3.23) 


X 


a(l — a){2a — l)Ao + 


In oririiEi/^n 


‘Sif {T - iaAo,s) = haiti,t 2 , s,b 2 , Xq) + 


1 /I - a \ .2 


2 V a 

with ha of ()3.2p . So, similarly to the proof of Lemma [2] one can write 

j^^nA+{l-a)\o{xi+X2) 


62 + (1 — Ci)^0 


F2{A) = Cn{X)- 


{xi - X2) 


x(^jg{T, + O ^ , (3.24) 


where Vn = U^-i/2yagn ((O) ^^2)) and A is defined in (13.191) . 
Repeating the argument of Lemma [21 we get 

F,(A) = n-’/=C„(.Y). ■■^"P{"^+(l-°Uo("l+"2)} ^ , 

[Xl - X2) 


/ ■")+5^™" ("■ V'’^) 


[— log n,log n]® 


X Q^U{T-iaAo,s)-A) 


where (•,•) denotes the scalar product in R^. The hrst integral is zero, because / is a 

symmetric with respect to ti and t 2 function and —— (O, ^^62 ) = —t;— (O, ^^ 52) = 1- 

ati ^ ^ ot2 
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Expanding the exponent into the Taylor series, we obtain 


F,(A) = „-V2c„(X). ^ 

[Xi - X2) 


X / fo, ^-^62^ g'^exp f-iaAo, ^-^62 ) dg. (3.25) 


It is easy to see that 




—2ixi i{xi — X 2 ) 0 

i{xi — X 2 ) 2 ix 2 0 

0 0 0 


Computing the integral in (I3.25h . we have (I3.17p . 

If Aq = 0, then a = As it was mentioned above, in this case there are two stationary 

points at the contour. The asymptotics of the integral (I3.2ip in the neighborhood of the 
second stationary point (i.e. (0,0,—!)) is computed by the same way as above. (|3.17p is 
transformed into ()3.20p . ■ 

The assertion (ii) follows from (|3.17ll and (I3.20p . 

Now we proceed to the proof of agreement between cases (i) and (ii) of Theorem [TJ 

Lemma 4 Let all conditions of Theorem[I\ are hold and Aq = 4 — 46^ — dn, dn —t 0. Then 
F 2 {A) satisfies the asymptotic relation 


F 2 {A) = y„exp{n(2 - 3bl - dn)/2 + Ao(xi + X2)/2}(I + o(l)). (3.26) 


where 




' -| /Q ^ 

ndrl , if dn > 0 ,ndn —t 00; 

< if nd'^ —)• const] 

^ C{-dn)~^/'^, if dn < 0, ndl ^ Qo. 


Proof. Consider the case dn > 0. Choose the same contour as in the proof of Lemma [2j 
Stationary points are also the same. 

Change of the variables t = ti + 12, a = ti — t 2 gives us 


F2{A) = CniX) 


jgAo(a;i+X2) 

2 (xi - X2) 


/ 




where f{T,a,s) = f{T,s). Set 

T* = -iXo] 

= (-1)1^4-462-Ag = (-irsr^ 

It is easy to see that 

fir*,a;,b2) = l-lbl-^dn] 


f (T*,cr*, 62 ) = - 


/ fei + <5„/4 {-l)^iXodn^^/4 ib2Xo/2 \ 

{-l)^iXodl/^/A dn!^ {-I)^b2dl/^ 12 
V ib2Xo/2 {-l)^b2di^^/2 bl + 1 / 


det f”{T*, a* b 2 ) = -dn/d] 
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{r*,a;,b2) = ^^^di/\5n-3y, 


d^f 

5^3V 

d^L * * , . 3 

^{r ,^„62 ) = -j. 

Let us choose Vn as a union of the products of the neighborhoods of r*, a\, 62 and 
T*, £72, 62 such that the radii of the neighborhoods corresponding to r and s are equal to 
\ogn/^/n, whereas the radius of the neighborhood corresponding to cr is equal to log nj^/n5n, 
if —>■ 00, and to log n/n^/^ otherwise. Similarly to the proof of Lemma [2] it can be proved 
that for (r, u, s) ^ Vn and sufficiently big n 


^f{T*,a*b2)-^f{T,a,s)>C 


log^ n 


n 


(3.27) 


Let n6n —^ 00. Then by the same way as before, with the only one distinction that the 
change of the variable u is fx —>■ + a*, we obtain 

ft(A) = n-“/2C„(X). + ^ 

2 {xi - X2) 

2 

^ loK 


1aM( x, 

X / ’ ’ drdads, (3.28) 

where is a quadratic form, defined by the matrix, which is obtained from /"(r*, cr*, 62) 
1/2 

by dividing by 5n of all numbers in the second line and the second column, i.e. 

= {T,6n^^‘^a,s)f"{T*,cr*,b2)iT,6n^^‘^a,sf. 


Therefore we have (|3.26p . 

Now let n5n —>■ 0. Then, changing the variables ci^ = d, we get 


+ 00 


pAQ[xi-\-x2) r r ~ j— 

F2m = Cn{X)- -- / dd / sin((xi - X 2 )Vd/ 2 )e-*(^l+^ 2 )r/ 2 gn/(r,A/S..)^^^^_ 

2(xi -X2) J J 

0 R2 

Let 1^ be a product of the ^^^-neighborhoods of 0 and 62 - Then we can shift the variable 


Vn 

T ^ T — i\o and in view of ()3.27l) restrict the integration domain by 


log^ n 
y/n 


X Vn 


F2{A) = Cn{xy 


pAo(3;i+a:2)/2 

2(xi - X2) 


y/n 


(1 + o{l))dd 


X [ sin((xi-X2)\/d/2)e-*(^i+^2)"/2gn/(r-*Ao,A/s,s)^^^^_ 


Vn 
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Expanding / and sin by the Taylor formula near (—zAq, 0, 62) and 0 respectively and changing 
the variables r a —>• s n~^l‘^s + b 2 , we have 


- 3 / 2 ^ , exp{n(2 - 3bl - (5„)/2 + Ao(xi + X2)/2} 
F2 (A) = n C„ (X)-2 {i.-X 2 )- 


log^ n 


j {l + o{l))da j (xi - (3.29) 


0 [—log n, log n]^ 

where is a quadratic form defined by the matrix 


2^ 


a^f 

1 

a^f 


2 drdcT^ 

drds 

1 a^f 

1 a\f 

1 a^f 

2 draY^ 

12^ 

2 

a\f 

1 aXf 

a^f 

ards 

2 dcr'^ds 



h ‘2 iAo/8 ib2y^ol‘i 

iAo/8 1/16 62/4 

y 1 , f .f j \ib2Xo/2 62/4 1 + 6| 

where all entries of the matrix P are units. Performing the Gaussian integration, we obtain 


■ exp{n(2 - 362 - <5n)/2 + Ao(xi + X2)/2} 

rv^!^ 4yn 

+00 ^ 

J Cdexp|-^|fid(l + o(l)) 

0 

that imply (13.261) . 

If n(5^ ^ const, then there is a certain third power polynomial instead of in the last 
exponent in (j3.29p . Thus, the asymptotics of F 2 {A) differs from (j3.26p only by multiplicative 
n-independent constant, with is absorbed by C in Yn- For negative 5n, if —)• 0, all the 

changes in equations appear in multiplication by factors which are equal to 1 + o(l), so (I3.26h 
is unchanged. If n5‘^ oo, then the combination of the argument of the case <5„ > 0 and of 
Lemma [3] causes some changes in (I3.28h and implies Yn = C(—Sn)~^''^ in (I3.26j) . ■ 


4 Proof of Theorem [2] 


As in the case of Theorem [H the proof of Theorem [2] is based on the application of the 
steepest descent method to the integral representation of F2m obtained in Section [2l For this 
end a “good contour” and stationary points of f 2 m, defined in (|2.8I) . have to be chosen. We 
start from the choice of stationary points. 

If p = n and hence bi = b 2 k = 0, Z > 1, the proper stationary points are 


where 


Set 



,j = l,2m, Ri = 0, 



T = diag{t*}2™ , 
b = (62) • • • ) b 2 mi ^ 4 ) ftg, . . . , b 2 m)- 


(4.1) 
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Since 


{A2m{T,R)f'^m{T,R))' 


T=T 

R=0 

6=0 


= (-1)™ + + = 1 or 2 

\2m) 

is nondegenerate for Aq G (—2,2), for sufficiently small b there exists the unique solution 

T = T{b), R = R{b) (4.2) 


of the equation 


A2miT,R)f^^{T,R) = 0 


(4.3) 


such that T(0) = T, R{0) = 0, and the solution continuously depends on b and Aq. When 
p —>■ oo, (12.lip and (|2.14p yield 6 —>• 0. Therefore, the solutions (14.2p are the required 
stationary points. 


Lemma 5 The solution (14.21) also has the following properties 

(V Tj^jiib) = Tj^j^{b), = 7fc2fc2(6) for all ji,j 2 G I+, ki,k 2 G I_, where 1+ = 

{j, t* = t*}, I_ = {j, t* = -t*}; 

( 2 ) {Ri)a 0 {b) = 0 , 1 = 2 , 2 m, aj^/3. 

Proof. Let 


TtT = diag{4(j)}, TtR = {ttR 2 , . . . , 7rR2m), TT G S 2 m, 


where ttRi is the such matrix that {'KRi)ai3 = {Ri)a.„j3.,,- Then Vvr G S 2 ra / 2 m( ttT, 7ri?) = 
f 2 m{T,R). So, it is sufficiently to proof the lemma for those stationary points, for which 


+* — — +* —_ +* — —_+* — f* 

A — ■ ■ ■ — ''k — '-fc+i — • • • — Am — ■ 


We are going to prove that there exists the solution of (14.31) that satisfies conditions 
(l)-(2) and r(0) = T, R{0) = 0. It is equivalent to existence of the solution of the system 


A 2 m{T,R)-^f 2 m{T,R) = 0, j = 1,2m; 


dt 


A2m{T, R) 


d 


(4.4) 


diRi)c 


-f 2 m{T, R) = 0, a € l2m,i, I = 2, 2m, 


where T and R satisfy (l)-(2), with respect to the variables A, t 2 m, {Ri)aa- Since the 
derivative of the l.h.s. of the system at T = T, i? = 0, 6 = 0 is nondegenerate, there exists 
the solution of it. In view of uniqueness of the solution of (j4.3h . the solution of (14.4p coincides 
with ()4.2I) . ■ 

The next step is the choice of the “contour” (in this case it is a d 2 m-dimensional man¬ 
ifold, d 2 m = (( 2 m) ~ — 1 + 2m)). For each variable we consider some contour and the 

required manifold M 2 m will be the product of these contours. Fix some variable. Order the 
corresponding components of the stationary points by increasing of the real part (if real parts 
are equal, order by increasing of the imaginary part). Then the contour is a polygonal chain 
that connect points by the order described above. Infinite segments of the polygonal chain 
are parallel to the real axis and directed from the first point to the left and from the last 
point to the right. 
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The Cauchy theorem and (j2.7p imply 


9 f zm 

j2m. -m |Ao ^ Xj | 

F^miA) = - / 92miT)e^^^-^^'^^dTdR, 

Al27n 

where 

f 2m 

g 2 m{T) = A(r) exp i - * ^ Xjtj 
^ i=i 

Moreover, 

F2m{A) = -( J 92m{T)e^^^-^^'^^dTdR + r{n,N) 

MiL 

where 

MiL = {C € M2m\ PCiloo < N}, 

|r(n, A)| < N ^ oo, 

llylloo = max|yA- 

3 _ 

Let e > 0 be an arbitrary positive number. Then, since {T(b),R{b)) —)■ (T, 0), for 

n^oo 

sufficiently big n and for every (T, R) E we have 


5R/L ( - AAo, ^R ] - ^fLiT, R) 


< e, 


where f^^{T,R) = f 2 m{T,R)\b=o- Also, f 2 m{T,R) ^ f^^{T,R), oo, {T,R) £ K for any 
compact set K. Hence, for sufficiently big n 

m 2 m{T, R) - ^fL{T, R)\<e, (T, R) E 

Consider the point {T^, R^) E such that = 5ftT, = 0. Then ^f 2 m{T^^ R^) > 
$R/o,,(f,0)-2e. Thus, 


max 

{T,R)eM^^ 


m2m{T,R)>mL{T,0)-2e 


max 5R/2m(r, 0) - 2e. 

^T=Ao/2 


Therefore, if 5R/ 2 m(T^, ii^) = max 3ft/2m(T, ii), then 


{T\R^) E {(T, R) E M^J 5R/2m(T, R) > mlm{T, 0) - 2s} 

C {(T, R) E M^J ^fLi^T - Uom > ^fLiT, 0) - 4e}. 

So, it is evident, that {T^,R^) —)• (T, 0), n —)• oo for certain T of the form (|4.ip . 

Let Vn{T{b), R{b)) be the neighborhood of the stationary point (T{b),R{b)), which 
contains the corresponding maximum point of 5R/2m with its ^^^-neighborhood, and 
diamlA(T(6), R{b)) —)• 0. It can be assumed that the union of these neighborhoods is invari¬ 
ant for map (T, R) —)• (vrT, ttR) for every vr E S 2 m- Then, by the same reasons as in the proof 
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of Theorem dl we can restrict the integration domain by the union of the neighborhoods Vn- 
Shifting the variables T —>• T + QT{b), i? —?■ i? + QR{b) in each neighborhood and expanding 
92 m by the Taylor formula, we get 


F2m{A) = 


■2m2-mexp 


f 2m N 

{ Ao E Xj I 

i=i 


x^e 


A{X) 

nf2^{T{b)Mb)) 


(1 + 0(1)) 

^ n-\-\/^D^g2m{T{h))F 

\a\<2m{m—l) 


Vr,{T{b),R{b)) 


where the summation is over all stationary points under consideration and 


VniT{b),Rib)) = Vn{T{b), R{b)) - 9(r(6), R{b)y, 

\{2m{m-l))/rp\\ ^ ^ u |2m(m-l)+l 
I 92m ^ '' — / ^ ' 3 ' 

3 

The number of terms of the Taylor series in (j4.5l] is the minimal number that allows us to 
obtain nonzero asymptotics. 

Fix some stationary point (T(5),i?(6)) and some multi-index a. Let /3 < a be a multi¬ 
index with = /3j2 for some ji / j 2 , ji ,^2 £ L+ or jij j 2 £ L_, where /? < a Vj I3j < aj. 
Then 


/ 

Vn{T{b)Mb)) 


^D“-^A(r)D^exp 



r 

T=T(b) 


+ D^-^A{T)D^ exp 





T=T{b) 


^nf2miT+^T{b),R+^R{b)) 


where a is the multi-index, which is obtained by swapping aj^ and aj^ in a. Hence, in the 
sum in (j4.5jl only the summands with |a| = 2m{m — 1) remain. 

Changing the variables T —)• R^ n we get 


j2m? m 


F2m(A) = n-'='2-/2c(2™)(X)- 


^ 2m N 

{Ao E Xj\ 
'' j=l 


A(X) 


(i + o(i))E 


,nf2miTib),R(b)) 


I (n—b ^ D^g2m{T{hW + ) 


Vn{T{b),R{b)) 


|Q|=2m(m—1) 


X exp -qf!y^{T{b), R{b))q^ M 1 + 


raiq) 


dq 


(4.6) 


where g is a vector which consists of all integration variables, dq = dTdR, Vn{T{b), R{b)) = 

^Vn{T{b),R{b)) and |r6(g)| <CEkiP- 

j 
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dMD implies that, as n ^ cx), 


{Rl^aoi 

(-^ 2 ) 00 : 


Tnib) 


0 ( 62 ), ^ = 3, 2m; 

_ &2 _ 

(o)r„ 2 „ 2 (o) 


+ 0 ( 62 ); 


aiai 


Tjj ( 0 ) 


( 2 m - x(j))rjj( 0 ) + xU) - 1 ^2 


Tni0r + Tj,{0f 


62 + o{b. 


where 


Therefore, 


xU) 


#/+, ifjG/+; 


Uf2m{T{b),R{b)) 


{-ir + xU){2m 


xU)) • 




+ 0 ( 62 ). 


Thus, the value of the ^f 2 m at the stationary points of the form (14.211 with #1+ = m is 
greater than that at the other stationary points of such a form for Aq G (— 2 , 2 )\{ 0 }. For 
Aq = 0 the values of the 3 ?/ 2 m at the stationary points are equal, because (j4.2p continuously 
depends on Aq- This yields that the sum in (I4.6P may be restricted to the sum only over the 
stationary points with #/+ = m (for Aq = 0 the other summands have the less order of n). 
We have 


•2m2-mexp 


F2™(A) = n-^2m/2-m{m-l)(j{2m)^j^y 


^ 2m N 

{Ao E Xj\ 

^ J =1 


A{X) 


-(1 + 0 ( 1 )) 


E ' 

#I+=m 


,nf2m(T{b),R{b)) 


I D^92m{T{b)Y 


|Q:|= 2 m(m— 1 ) 


X exp |-g/2m(r(6), R{h))q^ \ dq 

Consider the term with = {1, 2,... , m}. The Gaussian integration gives us 


Cri^ i 


2 .3^2_2^ A(xi, . . . , Xm)A(^Xm+l ; • • • ; X2m) 


f > 2m 

I 77177- , cy An V — > 

exp —(A2-2) + ^^a 
^ i=i 

X exp I - ~ f;(x,- - x™+,)|(l + 0(1)) = Cn™'(-l)™(-+i)/2 

^ i=i 


exp I - ^ 

X _ 

m 

iXi ®m+fc) 

iA=i 

where C is some real n-independent constant. 
On the other hand. 


exp<j ^(Ag -2) + ^'^Xj [>(1 + 0 ( 1 )), 


Aq 


2m 


i=i 


det 


S2miX) = 


^i7rpsc(^o)(^j ^m+k) _ g i'^PscAo}{^j ^m+k) 

{Xj Xm+fc) 


j,k=l 


(2i7rp5c(Ao))™'A(xi , . . . , X)7i)A(Xm+l) • • • J X2m') 
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where S 2 m is defined in (jl.7h and psc is defined in (II.4p . The determinant in the l.h.s. is the 
sumof exp{.:.,„(A„) over all collections {s,}. which consist of tnelements +1 and 


m elements —1, with certain coefficients. Since (see |19l Problem 7.3]) 


(- 1 ) 


m(m— 


n K' - Uk){vj - Vk) 


n (uj - Vk) 

j,k=l 


= det 


1 


'“'i j j,k=l 


the coefficient under exp < — iirpsci^o) ^ {xj — Xm+j) r is 

f j=i J 


det 


{Xm+j Xk) J j^k=l 




(2i7rpsc(Ao))™A(xi , . . . , XfYi ) , . . . , X2m) 


(2i7rpsc(Ao))™' D[ {xj Xjji^k) 
j,k=l 


The other coefficients can be computed by the same way. Therefore, 

2m 


F 2 m{A) = Cn^" exp <| ^ 


The assertion of the theorem follows. 


^ zrn N 

(Aq — 2) + — ^ Xj >,S 2 m(A)(l + o(l)). 
j=l ^ 


5 Proof of Theorem [3] 

In this section we consider the case p —)• oo. As in the proof of Lemma [3l the good contour is 
= ’At 2 = —aAo, s € M with the stationary point (—iaAo, —iaAo, ^^ 62 ) j where a satisfies 

(ICT) . 

Set /? = 2a — 1. Then (|3.18p transforms to 62 = /?(! + /?)(! — ■, and hence /3 = 

+ 0 ( 1 )). Substitution of Aq = ±2, a = and 62 = /3(1 + /3)(1 — /3)“^ in (I3.22p - (I3.23I) 

yields 


1 — a, 


a 


/3(2 + /3) 


f[-iaAo, - b 2 ]=-\ /3^(l-/3)(l + /3) 

±i/3(l - /3) 


-1 


1 — a, 
a 


4/3 


det /"( - iaAo, A—lb 2 ] = - , • (2 - (1 + /3)(1 - /3)(2 - /3)) = - 


/3i(l_/3)(l + /3)-i ±z/3(l-/3) 

/3(2 + /3) ±i/3(l-/3) 1; 

±i/3(l-/3) I + /32 

4/32(1+ 2/3-^2) 


(l + /3)^ 


(l + /3)^ 


In addition, 


^ T-iaAo, —^^2 ) = -2isign(Ao)(l - /3)^ 


2/3 1 

If (i) ^ -^ then Vn is chosen as a product of the neighborhoods of 0 , 0 , and ^-^b 2 such 

that the radius of the corresponding to ti and t 2 neighborhoods is but the radius of the 
corresponding to s neighborhood is . We have p.24p . 
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Change of variables T (n/3) s n + ^-^h 2 and repeating of the argument 

of the proof of Lemma [3] yield 

^2(2/ + n~‘^/^X) = Cpex.^{nA + n^/^Ao(xi + X2)/2}(1 + o(l)), 

where C is some absolute constant and A = f (—iaAo, ^^62). 

The assertion (i) follows. 

2 ^3 

Let now (ii) —>■ c. Chose Vn the same as in the case (i), but the radius of the 

neighborhoods corresponding to ti and t 2 is Then (I3.24p is also valid. 

The Cauchy theorem implies 


F2{2I + n-^/^X) = Cn{X) 


j^gnA+2n^/^{l—a)(xi+X2) 
n^/^{xi — X 2 ) 


Wr, 


+ 0 e 


,-Clog^n 


where the integration domain over s is not changed, but the ones over tj become {\z\ < 
n“^/^logn I argz = —tt/6 or arg^; = —57r/6}. 

Changing the variables T —?• s -A + ^-^62; we obtain 


F,_{21 + n-y^X) = n-'-f^CJX) ■ + 2^1^!- + x,)} ^ 

(rci - X 2 ) 


l+^f.,2 
X / e 2 ^ 


I' 


ds I (ti-t 2)exp<{ + ixjtj ) J^dT 


7 X 7 


j = l 


V2 


where 

7 = {arga: = —tt/G or arg^; = —Svr/G}. 

Change of the variables tj = tj + and using of the Cauchy theorem gives us 

F,i2I + n-^/^X) = C«=/= . »«P{"-^+(2"‘7l-°) + '/^)fa+"2)} ^ 

[Xi - X2) 

2 ^ . 

I 

3 


7x7 


(n - T2) exp <j - ^ ( -rf + i{xj + 2 c)Tj j |>dridT2, 
i=i 


where C is some absolute constant. Taking into account (ll.Sp and (j2.13p . we get 

^2(2/ + n~‘^^^X) = Cr?!'^ exp{nA + ( 2 n^/^(l — a) + V^)ixi + X2)} 

X A(xi + 2c, X 2 + 2c)(l + 0(1)). 


which completes the proof of the assertion (ii). 


6 Appendix 

6.1 Grassmann variables 

Consider the set of formal variables {V’j, which satisfy the anticommutation relations 

ipjipk + = 'il^jipk + '^k'll^j = = 0 - 
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This set generates a graded algebra A, which is called the Grassmann algebra. Taking into 
account that = ipj = 0, we have that all elements of A are polynomials of {tpj, 

We can also define functions of Grassmann variables. Let x be an element of A and / be 
any analytical function. By f{x) we mean the element of A obtained by substituting x ~ 
in the Taylor series of / near zq, where zq is a free term of x- Since x ~ is a polynomial 
of {'ipj, with zero free term, there exists I G N such that (x — zqY = 0, and hence the 

series terminates after a finite number of terms. 

The integral over the Grassmann variables is a linear functional, defined on the basis by 
the relations 



dipk = 0 , 




= 1 . 


A multiple integral is defined to be the repeated integral. Moreover “differentials” 
{d^lJj, anticommute with each other and with {'ijjj, Hence for a function / 


..., V’n) = «o + ^ + ... + n i’j 

i=i i=i 


we have by definition 

# 1=01 . 

The use of Grassmann variables for computing averages of determinants rests on the 
following identity, valid for any n x n matrix A: 



n ^ 

- \ n dtJjjdiljj 

j,k=l ^ j=l 


det A. 


( 6 . 1 ) 


The one more identity is the Hubbard-Stratonovich transformation 

. , ( 6 - 2 ) 

gj/i — ^ / ^ay{u+iv)+at{u-iv)-Av?-o?v'^ 

TT J 

which valid for any complex numbers y, t and any positive number a. The identities (j6.2h 
also hold when y, t are arbitrary even Grassmann variables (i.e. sums of the products of even 
number of Grassmann variables). For even Grassmann variables the formulas can be proved 
by Taylor-expanding and into the series and integrating each term. 

The properties explained so far suffice to obtain the integral representation for F 2 m, 
m = 1, whereas the general case m > 1 requires some additional preliminaries, pertaining to 
antisymmetric tensor products. Further details about antisymmetric tensor products may be 
found in 13Q1 Ghapter 8.4]. 



6.1.1 Grassmann variables and the exterior product 

The exterior product of vectors is well-known, as well as the exterior product of alternating 
multilinear forms (see m)- However, to prove Proposition [T] we need the exterior product of 
alternating operators. Define it as following. Let A be a linear operator on A'^C"' and B be 
a linear operator on A'’C"'. Then the exterior product A A H is the restriction of the linear 
operator Alto(A(8) B) on the A^+^C*^. Here Alt is the operator of the alternation, i.e., 

Alt(t) = ^ sgn7r/^(t), t € A^F, 

neSk 
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where Sk is the group of permutations of length k; sgnvr is the sign of permutation vr; is 
the canonical automorphism of which carries ui (8 )... <8) Ufc to u^(i) (S) ■ ■ ■ <S> Vj £ V; 
V is some finite-dimensional linear space. Note, that for A E End 18 the exterior product 
A A coincides with the well-known second exterior power of linear operator A. 

Fix some basis of Let A E EndA^C*^ and a,l 3 £ In^k, where In,k is defined 

in By we denote the corresponding entry of the matrix of A in the basis A 

... A , Oi E Iri^k\‘ 

To obtain the integral representation for F2m, we use the lemma: 

Lemma 6 Let A and B be linear operators on A'^C"' and A'’C” respectively. Then 


r 




0 !.yj 3 ^In,q j — 1 j — 1 



Proof. Let Sg^r be the set of such vr E Sg+r that satisfy inequalities 7r(l) < ... < 7r(g) and 
TT{q -|- 1) < ... < Tr{q + r). Then 


r 




OlA&In.q i = l 'y,S£l„,r j = l 


g+r 


Y X] vr sgn a, 



where 


OItt (Q^7r(l); • • • 5 ^7r(gH-r)); 

CX — (o^l; • • • ; ^q) ^ 

CX — (oig-i-i, • • • , CXq-\-r) G In^r- 


On the other hand, 


[A <S> A e^//) — (yA ® 






where Cq, = Ca^ A ... A Ca^, a E In,q- 
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Hence 


q\r\ ( \ 

Alt((H (g) B){epi A epu)) = - -yj- E E sgn aAapr^B^/s'J j A 

= ^q 'l'r)\ ^ ^ Ae„» 

O^^In^q-\-r TT^O'^Sq^r 


q\r\ 

{q + r)! 


E E sgn 7 rsgncrHa;; 3 ;Ha»/ 3 "e„/ A e^". 

CX.^In,q-\-r ^^^Sq^r 


Thus, 

- sgn7rsgncrHa;^^H„»;3», 

TT, (J G‘S'g, r 

which completes the proof of the lemma. ■ 

We also need some properties of the exterior product of the operators. 



Proposition 2 Let Aj G EndA'^^C"^, j = l,k, and B G EndC”. Then 

(i) Ai A A2 = A2 A Ai; 

(ii) (Ai A A2) A A3 = Ai A (A2 A As); 



k 

where q = Y] qj, B^^ = B A ... A B^. 

j=i ' "A ' 

q times 


Proof. Assertions (i) and (ii) follow from Lemma 1 and from Grassmann variables multi¬ 
plication’s anticommutativity and associativity. 

(iii) Consider the case k = 3 . 


Ai A A2 A A3 = Alt o((/ o Ai) (g) (Alt o(A 2 (g) As))) 

= Alt o(/ (g) Alt) o (Ai (g) A2 (g) A3) = Alt o(Ai (g) A2 (g) A3), 


where I is the identity operator. 

The general case follows from the induction, 
(iv) By definition, we have 


A B''i> A, 
i=i 


k 

Alto(g)H^ 9 ,A^., 

i=i 


29 



Consider (Alt ){vi ® ® Vq-)^ vi G C” 


(Alto5®®)(r;i 


' %■) = ~F X] sgn7r/^(Sr;i (g) ... ® Bvq/, 

Qj. 


T&Sn 


—j- sgnvri??;, 


(Ij 


'7r(l) 


I Bv. 




Tes„ 


^ ^ sgnvrS®® (/^(?;i (g) ... (g) -Ug^.)) = (5®® o Alt)(z;i (g) ... (g) Vg^). 


Tres,^ 

Therefore, Altoi?®'^i = B®'^^ o Alt, in particular, B^'^^ = |C" • Thus, 

k k y k 

f\ B^i^ A j = Alt o (g) B^^iAj = Alt oB®? o ( (g) Ag 
j=i i=i ^ i=i 

k 




The proof of the second formula is similar. 


6.2 The Harish-Chandra/Itsykson—Zuber formula 

For computing the integral over the unitary group, the following Harish-Chandra/Itsykson- 
Zuber formula is used 


Proposition 3 Let A be a normal n x n matrix with distinct eigenvalues and B = 

diag{bj}'j^i, bj G M. Then 


J exp{ztrAU*BU}dUn{U) 



det{exp{zajbk)}'J,,^i 
zX^-^y^A{A)A{B) ’ 


where z is come constant, A(A) = A({ag}”^g), A is defined in (I 2 . 1 h . Moreover, for any 
symmetric domain H and any symmetric function f{B) of {bj}^^-^ 


exp{ztrAU*BU}A‘^{B)f{B)dUn{U)dB 


Un ^ 


^ •- „-(n2-n)/2 

I ■ ^ 


i=i 


A(A) 


J exp |2;^ag6g |A(B)/(B)dB, ( 6 . 3 ) 
n 


n 

where dB = n dbj. 
i=i 

For the proof see, e.g., [HI Appendix 5 ]. 


Remark 2 Notice, that (16.3p is also valid if A has equal eigenvalues. 


Indeed, if Ogj = aj^, ji 7^ 32-, the integrand in ( 16 . 31 ) changes sign when bj.^ and hj^ are swapped. 
Thus, the ratio of the integral and Ug^ — Ug^ is well defined. 
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